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Abstract 



Nonrelativistic bound states have two low energy scales, a soft scale us 
of order mv and an ultrasoft scale \x\j of order mv 2 . In two-stage running, 
the soft and ultrasoft scales are lowered from m to mv, and then the ultrasoft 
scale is lowered from mv to mv 2 . In one-stage running, the two scales are low- 
ered in a correlated way using a subtraction velocity. We compare these two 
methods of summing logarithms and show that only the correlated running 
in velocity space of the one-stage method correctly reproduces the logarithms 
in nonrelativistic bound states in QED. The argument for one-stage running 
is general, and should apply to any system with correlated scales. 
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Typeset using REVTf^X 



Nonrelativistic bound states have two important scales, a soft scale [is of order mv and 
an ultrasoft scale [iu of order mv 2 , which are of order the typical momentum and energy of 
the bound state In a Coulombic bound state such as Hydrogen, the typical velocity 

v is of order the coupling constant a. Logarithms of v or equivalently a can be summed 
using renormalization group equations (RGE) This is accomplished by constructing 

an effective theory below the scale m that contains soft and ultrasoft modes PHl^f, and then 
scaling the coefficients in the effective Lagrangian to low energies using the RGE. In a prop- 
erly constructed effective theory, all the large logarithms should arise from renormalization 
group running, and there should be no large logarithms in the matching conditions. In QCD 
the summation of logarithms is important because a s In v ~ 1 and the RGE systematically 
accounts for the running of a s . In QED, a In a is small but higher order corrections are 
necessary to confront precision experiments. 

In two-stage running, one first lowers the subtraction scale \x from m to mv. At this point 
the soft modes, with energy- momentum of order mv, are integrated out of the theory. The 
ultrasoft modes are then lowered from mv to mv 2 . This is the definition of two-stage running 
that we will use in this article. [| In one-stage running [[|, one introduces a subtraction 
velocity v, with [i$ = mv and \ijj — vnv 2 . The subtraction velocity v is then lowered from 
v — 1 to v — v, so that lis — mv an d [lu — mv<2 & t the end of the renormalization group 
evolution. The difference between the two methods can be seen in Fig. [1], where the path in 
the Hs — ^u plane is shown. In contrasting the two methods we will consider NRQED, where 
the analysis is simplified by the fact that the coupling constant a does not run below the 
electron mass. Both the two-stage and one-stage methods correctly reproduce observables 
with a single In a, such as the a 5 In a Lamb shifts in Hydrogen and positronium |TB|fr7| ,[5[] . In 
Ref. H one-stage running was shown to correctly reproduce the a 7 In 2 a hyperfine splittings 
for Hydrogen, positronium and muonium [|l8| -|2"0|j , the a 3 In 2 a ortho and para-positronium 
widths |[l|j2gj2|, and the a 8 In 3 a Hydrogen Lamb shift |]|-0 This 

is a highly non-trivial 

check of the one-stage method, and is the first time higher-order logarithms in QED bound 
states were computed using the renormalization group. In previous computations, higher 



order logarithms were obtained by direct computation of multi-loop diagrams |18|-t22|. In 



this article, we explain why the two-stage method fails to reproduce the In 2 a and In 3 a 
terms, and why the one-stage method is essential for correctly summing all the logarithms 



See, for example, Sec. 1.4 of Ref. ||15|| . The running in Ref. corresponds to the m to mv 
running of the two-stage method. 



2 Note that the a 8 In 3 a term has also been calculated numerically in Refs. p3|-p5| and 
there is general agreement on the value of the 4-loop diagrams calculated analytically by 

0. 



Karshenboim in Ref. 
another graph. 



However, Refs. ||23|J25|1 have an additional contribution from 
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FIG. 1. Paths in the fijj — us plane for one-stage and two stage running 



using the renormalization group equations. 

A loop graph computed in the effective theory below m typically contains logarithms 
of the form lnp//ig, In Ej flu and log y/mE/jjs, where E and p are the typical energy and 
momentum in the bound state. We will define the soft and ultrasoft anomalous dimensions 
7 5 and 7^ by differentiating with respect to fis and pu respectively. Introducing both ps 
and nu allows us to simultaneously discuss the one and two-stage running. In the effective 
theory the potential is the Wilson coefficient of a four-fermion operator and has an expansion 
in v. The coefficients {Ui} in the potential have anomalous dimensions: 



7f, 



(1) 



dp s ' r " d/iu 

where j s and 7^ are functions of a and potential coefficients {Uj}. 

In two-stage running, one sets \x = p,$ = Hu and lowers \x from m to mv. Thus, one runs 
the potentials between m and mv using the anomalous dimension V s ^ 

da 



/'■ 



d/x 



r? 



7f + 7f 



mv < p, < m . 



(2) 



One then integrates out the soft modes at p — mv, and runs to mv 2 with the ultrasoft 
anomalous dimension 



3 In pNRQED, Tf is called the soft anomalous dimension because it governs the running 
between the hard scale m and the soft scale mv. In the one-stage approach, it is 7,- which 
is called the soft anomalous dimension. This is merely a difference in terminology with no 
physical consequences. 
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^— - = ji mv 2 < < mv . (3) 

In one-stage running the subtraction scales fxs and \iu are different ,Q/i,s = mv and nu = mv 2 . 
With the velocity renormalization group one runs from v = 1 to v = v with the anomalous 
dimension 

^ = 7 f + 2 T f v<v<l. (4) 

As an example, consider the running of the 5-function contact potential for Hydrogen 
in the limit m proton — > oo. The coefficient of this potential will be defined as U 2 , using the 
notation of Ref. ||. U 2 plays a key role in NRQED since it is the lowest order coefficient 
that runs. With one-stage running U 2 has a velocity renormalization group equation of the 
form 

v^ = 2a 2 \ u + a\ s Ul + ... . (5) 
dv z 

Here X u = — 4/(3m 2 ), A 5 = m 2 /iT, and the matching value U^iy = 1) = 7ra/(2m 2 ), where 
m is the electron mass. The coefficient 11% is of order a, so the first term in Eq. (|5]) is 
the leading order (LO) anomalous dimension, and the second term belongs to the next-to- 
leading order (NLO) anomalous dimension. For Hydrogen, the first term is purely ultrasoft 
and the second term is purely soft, so 
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u - a 2 X u , 7 5 = a\ s U 2 . (6) 



In Eq. (|5|), we have only retained the NLO term relevant for the a 8 In 3 a Lamb shift in 
Hydrogen. The full equation is given in Ref. § and gives additional contributions to a 7 In a 
and a 6 In a, but the simplified version in Eq. @ will suffice for this paper. 

In our example, integrating the renormalization group equation using the LO anomalous 
dimension and one-stage running gives 

U 2 (y = v) = U 2 (l) + 2a 2 X u \nv, (7) 

which is the leading order series a n+1 ln n v for U 2 {v) and determines the a 5 In a Lamb shift 
correction. There is only a single-log term since the LO anomalous dimension is a constant. 
To next-to-leading log order only the leading log value of U 2 is needed on the RHS of Eq. 
Substituting Eq. (|7|) for U 2 into Eq. (|5|) and reintegrating gives 



4 In dimensional regularization where d = 4 — 2e, factors of /x| and appear along with 
interaction vertices involving soft and ultrasoft modes respectively 
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U 2 (y = v) = U 2 {1) + 2a 2 X u \nv 



+aX ti U 2 (lj 2 \nv + 2a 3 X u X s U 2 (l)\n z v + (\ u )~ \»\n 6 v 



i v 



(8) 



which includes the LO series a n+1 ln n v and the next-to-leading order series a n+2 ln n v. The 
highest logarithm in the NLO series is ln 3 t>, since the NLO anomalous dimension has a U 2 
term that contains In 2 v terms when the LO running for U 2 is used. The In 3 v term gives an 
a 8 In 3 a Lamb shift for Hydrogen in agreement with Karshenboim ||18||. 



Next we redo this example using two-stage running. In the first step, one integrates 

dU 2 



dfj, 



a z X u + aX a U. 



2 J 



(9) 



to give 



U 2 (fi = mv) = U 2 (l) + a 2 X u lnv 

+a\ s U 2 (l) 2 \nv + a 3 X u X s U 2 (l) In 2 v + ^a 5 (X u ) 2 X s In 3 

where U 2 {1) = U 2 {ji = m) . In the second step one integrates 

dU 2 



d/i 



a 2 X u . 



(10) 



(11) 



Since the X s term comes solely from the soft scale it does not appear in the RGE for the 
second step. Integrating Eq. (|TTf ) using Eq. fllCf ) as the starting point gives the final result 



U 2 (fi = mn' 



U 2 (l) + 2a 2 X u \nv 



+aX s U 2 (l) 2 In v + a 3 X u X s U 2 (l) In 2 v + -a 5 (X u f X s In 3 



v . 



;i2) 



Comparing Eqs. (§), and Eqs. (|1"2"D, we see that the final expressions agree in the single- log 
term, but disagree in the higher-order logarithms. 

One can study the difference between one- and two-stage running in the general case, by 
integrating the renormalization group equations perturbatively. The single-log terms (i.e. 
proportional to lnt> ) are obtained by integrating the renormalization group equations using 
the values for the anomalous dimensions at the matching scale, i.e. by substituting the tree- 
level values for {Ut} in 7 ({Ui}), which will be denoted by 7(0). In the two-stage approach, 
one finds that 



Ui(n = mv 2 ) = rf(0)lnt; + 7f(0)lnt; = 7 f(0) + 7l u (0) Inv + 7^(0) lnv 



Ut 



(13) 



where the first term is the contribution between m and mv, and the second term is the 
contribution between mv and mv 2 . In the one-stage approach, one finds 



Ui{y = v)= 7 f(0) + 2 7 f(0) 



lnt>. 



(14) 
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Both methods give the same answer for the lnt> term, since 

If (0) + 7 f (0) = [ 7 f (0) + if (0)] + 7 f (0) = [ 7 f (0) + 2 7 f (0)] . (15) 

This is to be expected, since the anomalous dimensions in both methods should reproduce 
the single logarithms in the graphs from which they were calculated. This is true of our 
Hydrogen example, cf. Eq. (|8p and Eq. (|T2"D. We have also checked that this is true for the 
a 5 In a Lamb shift for positronium in which case the LO anomalous dimension has both an 
ultrasoft and soft contribution. 

The non-trivial terms are the higher order logarithms, which are generated when one 
includes the running potentials in the expression for the anomalous dimensions / y s,u . The 
anomalous dimensions depend on t = Inv through their dependence on L^. Consider ex- 
panding 7 in powers of t, 

where the derivatives with respect to t can be evaluated using the chain rule, 

13 [ ) 13 1 J \dU k dt ) t=Q + 2 [dUM dt dt + dU k dt 2 ) t=0 + • • ' 

\dU k J t=0 2 \dU k dUi dU k dUi J t=Q 

For Hydrogen, the In 2 v and In 3 v terms are generated from cross-terms of the form 
jU QrySjQjj^ an( j 0y u ^ d 2n f S /dU% that involve both soft and ultrasoft anomalous dimen- 
sions. In the one-stage approach, the renormalization group equations are integrated with 
the anomalous dimension in Eq. @ between v = 1 and v = v, so the cross-terms have the 
structure 



(2 7 ^)gln 2 ,, (2 7 ^) 2 ^ln 3 ,. (18) 



In the two-stage approach, one integrates Eq. (^j from m to mv, after which 7 s = 0, so the 
corresponding terms have the structure 



(l U )^T^v, (7T^m^, (19) 



d l S 1 2 / u\2 #V 



which differ by a factor of two and four, respectively, from the one-stage running values.^ 



5 Recall that there are terms in the NLO anomalous dimension that are not shown in Eq. ([5]) 
which contribute to terms involving ln a in the energy. Including these additional terms 
makes the first term in Eq. ( |18|) and in Eq. ([19|) more complicated, however the one- and 
two-stage In 2 v results still differ by an overall factor of 2. 
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FIG. 2. A four-loop graph that contributes to the a 8 In 3 a Lamb shift in Hydrogen. 

In general, one-stage and two-stage running give different values for the final answer. 
Define the anomalous dimension vector 7 = (7 C7 ,7 5 ). The running is given by integrating 
the anomalous dimension vector along the path shown in Fig. |I[ The integral is path 
dependent if V x 7 ^ 0, i.e. 

^f. + *£. (20) 
In our Hydrogen example, 7^ = a 2 X u and 7 s = aX £/f , so that 

^pL = 2a X s U 2 fiu^- = 2a 3 X s X u U 2 , 

dfjLu dfjLu 

d-y u 

tof- = 0, (21) 

and the integral is path-dependent. 

In Ref. ||, one-stage running was shown to correctly give the a 8 In 3 a Lamb shift and 
a 7 In 2 a hyperfine splitting for Hydrogen, muonium and positronium, as well as the a 3 In 2 a 
correction to the decay width for positronium. Since the answer differs from two-stage 
running, this implies that two-stage running does not sum all the logarithms of v. The 
reason can be seen by considering a diagram with mixed soft and ultrasoft divergences, such 
as the four-loop graph in Fig. [2|. If we calculate logarithms of a directly by calculating matrix 
elements (using vertices at the hard scale), then this graph contributes to the a 8 In 3 a Lamb 
shift. Calculating either ultrasoft subgraph in Fig. [| gives a term of the form ln(k°/nu), 
whereas the two-loop potential subgraph gives a term of the form \n(Vmk° / us), where k° is 
the time-component of a loop momentum. In the RGE, to make sure that both logarithms 
are small it is essential that \xy = [I'jg/m, so that the soft and ultrasoft scales are correlated. 
If one attempts to set [is = [iv — then one cannot simultaneously minimize both soft 
and ultrasoft logarithms inside this four-loop graph. In the two-stage trajectory shown in 
Fig. [I], the logarithms \n(k°/fiu) and ln(vmF//is) are made small at the endpoint of the 
path. However, both logarithms are not small everywhere along the integration path, so that 
the logarithms are not minimized inside subgraphs of the full graph. As a result, two-stage 
running gives the correct single log, but not the higher order logarithms. In the one-stage 
approach the logarithms are small along the entire integration path and the higher order 
logarithms are correctly reproduced. 
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Although two-stage running seems natural within the pNRQED framework fT0||i6| . |T7|j , a 
single stage running can be implemented as well. The discussion above suggests that even 
though matching from QED to pNRQED can be done in a two-stage procedure through 
NRQED, the renormalization group improvement should be done in one stage. However, it 
may well be that a more sophisticated version of two-stage running does properly take into 
account the fact that the soft and ultrasoft scales are correlated. 

The difference between one- and two-stage running is a generic feature of systems with 
correlated scales. In nonrelativistic bound states, the soft and ultrasoft scales are determined 
in terms of the velocity v, with /x s = mv and /jljj = mv 2 , and one needs to use one-stage 
running in velocity space. We have illustrated this point using QED, but the arguments are 
obviously also valid for weak coupling nonrelativistic QCD bound states. Correlated scales 
also occur in the problem of Sudakov logarithms near kinematic endpoints (for effective 
theory approaches see Refs. f26|). In B — ► X s ^ decay, there are two correlated scales, the 
collinear scale vn^X — x and the soft scale m&(l —x) (here x = 2E^/m,b) that are important 
as x — > 1. To sum subleading Sudakov logarithms using the renormalization group one needs 
to run in x (or equivalently moment) space. One-stage running might also be applicable to 
finite temperature field theory, where one has the correlated scales T, gT and g 2 T. This will 
be discussed elsewhere. 
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